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Abstract 

We investigate the connection between local and global dynamics of 
two Af-degree of freedom Hamiltonian systems with different origins de- 
scribing one-dimensional nonlinear lattices: The Fermi-Pasta-Ulam (FP- 
U) model and a discretized version of the nonlinear Schrodinger equation 
related to the Bose-Einstein Condensation (BEC). We study solutions 
starting in the immediate vicinity of simple periodic orbits (SPOs) rep- 
resenting in-phase (IPM) and out-of-phase motion (OPM), which are 
known in closed form and whose linear stability can be analyzed exactly. 
Our results verify that as the energy E increases for fixed A'^, beyond 
the destabilization threshold of these orbits, all positive Lyapunov expo- 
nents Li, i — 1, . . . , TV — 1, exhibit a transition between two power laws, 
Li oc E'^'' , Bfe > 0, k = 1,2, occurring at the same value of E. The 
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destabilization energy Ec per particle goes to zero as A'^ ^ oo following 
a simple power-law, Ec/N oc N'"", with a being 1 or 2 for the cases 
we studied. However, using the SALI, a very efficient indicator we have 
recently introduced for distinguishing order from chaos, we find that the 
two Hamiltonians have very different dynamics near their stable SPOs: 
For example, in the case of the FPU system, as the energy increases for 
fixed N , the islands of stability around the OPM decrease in size, the orbit 
destabilizes through period-doubling bifurcation and its eigenvalues move 
steadily away from — f, while for the BEC model the OPM has islands 
around it which grow in size before it bifurcates through symmetry break- 
ing, while its real eigenvalues return to +1 at very high energies. Further- 
more, the IPM orbit of the BEC Hamiltonian never destabilizes, having 
finite— sized islands around it, even for very high N and E. Still, when 
calculating Lyapunov spectra, we find for the OPMs of both Hamiltonians 
that the Lyapunov exponents decrease following an exponential law and 
yield extensive Kolmogorov-Sinai entropies per particle Kks/N cx const., 
in the thermodynamic limit of fixed energy density E/N with E and 
arbitrarily large. 

Keywords: Hamiltonian systems, Simple Periodic Orbits, regular and chaotic 
behavior, Lyapunov spectra, Kolmogorov entropy, SALI method. 



1 Introduction 

Chaotic behavior in Hamiltonian systems with many degrees of freedom has been 
the subject of intense investigation in the last fifty years, see e.g . [Lichtenberg fc 
Lieberman, 1991; MacKay & Meiss, 1987; Wiggins, 1988] and |Sim6 ed.. 1999j . 



By degrees of freedom (dof) we are referring to the number of canonically con- 
jugate pairs of positions and momentum variables, qk and pk respectively, with 
fc = 1, 2, . . . , iV. The relevance of these systems to problems of practical concern 
cannot be ov eremphasized. Thei r applications range from the stability of the 
solar system Contopoulos. 2002ll and the containment of char ged particles in 



high intensity magnetic fields [Lichtenberg fc L ieberman, 199l| to the blow-up 
of hadron beams in high energy accelerators [Scandale fc Turchetti. eds. 199l| 
and the underst anding of the properties of simple molecules and hydrogen- 
bonded systems [Bountis ed.. 1992t IProsmiti fc Farantos. 1995j . 

One of the most fundamental areas in which the dynamics of multi-degree of 
freedom Hamiltonian systems has played (and continues to play) a crucial role 
is the study of transport phenomena in one-dimensional (ID) lattices and the 
role of chaos in providing a link between determ i nistic and s tatistical behavior 
" Chirikov. 19791 iLichtenberg fc Lieberman. 199lt iFord. 1992j . In this context, 



a lattice of N dof is expected, in the thermodynamic limit {N ^ oo at fixed 
energy density E/N) to exhibit chaotic behavior for almost all initial configu- 
rations, satisfying at least the property of ergodicity. This would allow the use 
of probability densities, leading from the computation of orbits to the study of 
statistical quantities like ensemble averages and transport coefficients. 

Chaotic regions, where nearby solutions diverge exponentially from each 
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other, provide an excellent "stage" on which such a desired transition from clas- 
sical to statistical mechanics can occur. However, the presence of significantly 
sized islands (or tori) of quasiperiodic motion, in which the dynamics is "stable" 
for long times, preclude the success of this scenario and make any attempt at 
a globally valid statistical description seriously questionable. These tori occur 
e.g. around simple periodic orbits which are stable under small perturbations 
and, if their size does not shrink to zero as N or E increases, their presence can 
attribute global consequences to a truly local phenomenon. 

That such important periodic orbits do exist in Hamiltonian lattices, even 
in the N —^ oo limit, could not have been more dramati cally manifested than in 
the remarkable discovery of discrete breathers (see e.g. Flach fc Willis. 1998j '). 



which by now have been observed in a great many experimental situations 
[Eisenberg et ai, 1998; Schwarz et al., 1999; Fleischer et ai, 2003; Sato et 
ai, 2003]. Discrete breathers are precisely one such kind of stable periodic or- 
bits, which also happen to be localized in space, thus representing a very serious 
limitation to energy transport in nonlinear lattices. 

Then, there were, of course, the f amous numerical e xperiments of Fermi, 
Pasta and Ulam of the middle 1950's Fermi et al. 1955| . which demonstrated 



the existence of recurrences that prevent energy equipartition among the modes 
of certain ID lattices, containing nonlinear interactions between nearest neigh- 
bors. These finite, so-called FPU Hamiltonian systems were later shown to 
exhibit a transitio n to "global" chaos, at hi gh enough energies where major 
resonances overlap Izrailev fc Chirikov. 19661 . Before that transition, however, 



an energy threshold to a "weak" form of chaos was later discovered that relies 
on the interaction of the first few lowest frequency modes and, at least for the 
FPU system, does appear to ensure equipartition among all modes [De Luca et 
ai, 1995; Dc Luca fc Lichtenberg, 2002]. Interestingly enough, very recently, 
this transition to "weak" chaos was shown to be closely related to the desta- 
bilizati on of one of the lo west frequency nonlinear normal mode of this FPU 
system Flach et ai. 2005j . Thus, today, 50 years after its famous discovery, the 



Fermi-Pasta-Ulam problem and its transition from r ecurrences to true statisti - 
cal behavior is still a subject of ongoing investigation [Herman fc Izrailev. 2004j . 

In this paper, we have sought to approach the problem of global chaos in 
Hamiltonian systems, by considering two paradigms of N dof, ID nonlinear 
lattices, with very different origins. 

One is the famous FPU lattice mentioned above, with quadratic and quartic 
nearest neighbor interactions, described by the Hamiltonian 



1 ^ ^ /I 1 \ 



where Xj is the displacement of the jth particle from its equilibrium position, Xj 
is the corresponding canonically conjugate momentum of Xj, /3 is a positive real 
constant and E is the value of the Hamiltonian representing the total energy of 
the system. 
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The other one is obtained by a discretization of a partial differential equation 
(PDE) o f the nonlinear Sch rodinger type referred to as the Gross-Pitaevskii 
equation [Dalfovo et al. 199 9], which in dimensionless form reads 

ih^^ = --—^^ + V{x)^>{x,t)+g\^{x,t)\^^>{x,t), i^ = -l (2) 

where h is the Planck constant, g is a positive constant (repulsive interactions 
between atoms in the condensate) and V{x) is an external potential. 

Equ ation P]l is related t o the phenomenon of Bose-Einstein Condensation 
(BEC) [Ketterle et al.. 199^ . Here we consider the simple case V{x) = 0, fi = 1 
and discretize the x-dependence of the complex variable ^{x,t) = j{t) in 
(|2Jl, approximating the second order derivative by "^^x — ■ Set- 

ting then = Qjit) + i ■ pj{t), = -1, j 1, 2, . . . , iV and = 

q]{t) + Pj{t), one immediately obtains from the above PDE (O a set of or- 
dinary differential equations (ODEs) for the canonically conjugate variables, 
Pi and qi, described by the BEC Hamiltonian jTrombettoni fc Smcrzi. 20()li: 
ISmerzi fc Trombettoni. 2003j 

N N N 

j=i i=i j=i 

where 7 > and e ~ 1 are constant parameters, g = ^ > with Sx ~ 1 and E 
is the total energy of the system. 

In Sec. 12 wc study these Hamiltonians, focusing on some simple periodic or- 
bits (SPOs), which arc known in closed form and whose local (linear) stability 
analysis can be carried out to arbitrary accuracy. By SPOs, we refer here to 
periodic solutions where all variables return to their initial state after only one 
maximum and one minimum in their oscillation, i.e. all characteristic frequen- 
cies have unit ratios. In particular, we examine first their bifurcation properties 
to determine whether they remain stable for arbitrarily large E and TV, having 
perhaps finitely sized islands of regular motion around them. This was found to 
be true only for the so-called in-phase-mode (IPM) of the BEC Hamiltonian 

SPOs corresponding to out-of-phase motion (0PM) of either the FPU ^ 
or the BEC system destabilize at energy densities jf- oc N~", with a = 1 or 
2 (for the SPOs we studied), as ^ cxd. The s a.me result was also ob tained for 
what we call the OHS mode of the FPU system Oovama et al.. 1969j , where all 
even indexed particles are stationary and all others execute out-of-phase oscil- 
lation, under fixed or periodic boundary conditions. All these are in agreement 
with detailed analytical results obtained for families of SPOs of the same F PU 
system under periodic boundary conditions (see e.g. Poggi fc Ruffo. 19*9^ 



Then, in Sec. |31 we vary the values of E and N and study the behavior of 
the Lyapunov exponents of the OPMs of the FPU and BEC Hamiltonians. We 
find that, as the eigenvalues of the monodromy matrix exit the unit circle on 
the real axis, at energies Q < Ec = Ei < E2 < . ■ for fixed N , all positive 
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Lyapunov exponents Li, i = l,...,7V— 1, increase following two distinct power 
laws, Lj (X E^'' , Bk > 0, k = 1,2, wit h the i?^'s as reported in the literature, 
Rechester et ai. 19791: iBenettin. 1984 ^^'^ Livi et ai. 1986j . Furthermore, as 



the energy E grows at fixed N, the real eigenvalues of the 0PM orbit of FPU 
continue to move away from —1, unlike the 0PM of the BEG Hamiltonian, 
where for very large E all these eigenvalues tend to return to +1. Interestingly 
enough, the IPM of the BEG Hamiltonian remains stable for all the energies 
and number of dof we studied! 

In Sec. 01 we turn to the question of the "size" of islands around stable 
SPOs and use the Smaller Alignment Index (SALI) , introduced in earlier pap ers 
" Skokos. 200itlSkokos et ai. 2003at Iskokos et al. 2003bHSkokos et ai. 200^ to 



distinguish between regular and chaotic trajectories in our two Hamiltonians. 
First, we verify again in these multi-degree of freedom systems the validity of the 
SALI dependence on the two largest Lyapunov exponents Li and L2 in the case 
of chaotic motion, to which it owes its effectiveness and predictive power. Then, 
computing the SALI, at points further and further away from stable SPOs, we 
determine approximately the "magnitude" of these islands and find that it van- 
ishes (as expected) at the points where the corresponding OPMs destabilize. In 
fact, for the 0PM of the FPU system the size of the islands decreases mono- 
tonically before destabilization while for the BEG orbit the opposite happens! 
Even more remarkably, for the IPM of the BEG Hamiltonian ||2Jl, not only does 
the "size" of the island not vanish, it even grows with increasing energy and 
remains of considerable magnitude for all the values of E and N we considered. 

Finally, in SecEl using as initial conditions the unstable SPOs, we compute 
the Lyapunov spectra of the FPU and BEG systems in the so-called thermo- 
dynamic limit, i.e. as the energy E and the number of dof N grow indefinitely, 
with energy density ^ fixed. First, we find that Lyapunov exponents, fall on 
smooth curves of the form Li w Lie~"*/^, for both systems. Then, computing 
the Kol mogorov-Sinai entropy hr^c; , as the sum of the positive Lyapunov ex- 
ponents [Pesin. 1976HHilborn. 199l |. in the thermodynamic limit, we find, for 
both Hamiltonians, that Hks is an extensive quantity as it grows linearly with 
N, demonstrating that in their chaotic regions the FPU and BEG Hamiltonians 
behave as ergodic systems of statistical mechanics. 
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2 Simple Periodic Orbits (SPOs) and Local Sta- 
bility Analysis 

2.1 The FPU model 

2.1.1 Analytical expressions of the OHS mode 

We consider a ID la ttice of N particles with nearest neighbor interactions given 
by the Hamiltonian |Fermi et al.. IQSSj 

1 ^ ^ /I 1 \ 

^ = o E + E 9(^^+1 - ^^o' + i/'C^.+i - ^jf (4) 



j = l j=0 ^ 

where Xj is the displacement of the jth particle from its equilibrium position, Xj 
is the corresponding canonically conjugate momentum of Xj and /3 is a positive 
real constant. 

Imposing fixed boundary conditions 

Xoit) =XN+l{t) =0,'^t (5) 
one finds a simple periodic orbit first studied by Oovama et al. 1969j , taking N 



odd, which we shall call the OHS mode (using Ooyama's, Hirooka's and Saito's 
initials) 

iV - 1 

X2j{t) ^ 0, X2j-l{t) = -X2j + l{t) = X{t), j = 1, . . . , . (6) 

Here, we shall examine analytically the stability properties of this mode and 
determine the energy range < £' < Ec{N) over which it is linearly stable. 
The equations of motion associated with Hamiltonian 10} are 

Xj{t) = Xj+i-2xj+Xj-i+ p(^{xj+i - Xjf ~ {xj ~ Xj^if^ , j = l,...,N (7) 

whence, using the boundary condition jsj and the expressions (|SJ for every 
J = 1, 3, . . . , iV — 2, iV, we arrive at a single equation 

S:{t) = -2x{t) ~ 2Px^{t) (8) 

describing the anharmonic oscillations of all odd particles of the initial lattice. 
The solution of is, of cours e, well-known in terms of Jacobi elliptic functions 



Abramowitz fc Steeun. 1965j and can be written as 

x{t) ^Ccn{Xt,K^) (9) 

where 

and is the modulus of the cn elliptic function. The energy per particle of this 
SPO is then given by 
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2.1.2 Stability analysis of the OHS mode 

Setting Xj = xj + yj in {Tj) and keeping up to linear terms in yj we get the 
corresponding variational equations for the OHS mode (|SJ) 

y, = (1 + 3(3x^)iyj-i - 2y, + y,+i), j^l,...,N (12) 

where yo = yN+i = 0. 

Using the standard method of diagonalization of linear algebra, we can sep- 
arate these variational equations to N uncoupled independent Lame equations 
[Abramowitz & Stcgun, 1965] 

z^{t)+4{l + 3Px^)sm'(^-^^^^y,{t) = 0, j = l,...,N (13) 

where the Zj variations arc simple linear combinations of i/j's. Using and 
changing variables to u = At, Eq. (|13() takes the form 

z'Jiu) + 2{l+4K^-6Khii^{u,K^))sm'(^ ^^^^^ = 0, J = l,---,N (14) 

where we have used the identity cn^(u, k^) = 1 — sn^(w, k^) and primes denote 
differentiation with respect to u. 

Equation p4() is an example of Hill's equation [Copson, 1935; Magnus & 
Winkler, 1966] 

z"{u) + Q{u)z{u) = (15) 

where Q{u) is a T-periodic function (Q(u) = Q{u + T)) with T = 2JC and 
K. = /C(k^) is the elliptic integral of the first kind. 

According to Floquet theory [Magnus fc Winkler. 196(1 the solutions of Eq. 
p5|l are bounded (or unbounded) depending on whether the Floquet exponent 
a, given by 

cos(2a/C(K2)^ = 1 - 2sin2(A:(K2)^^ dct(D(0)) (16) 

is real (or imaginary). The matrix D(a) is called Hill's matrix and in our case 
its entries are given by 

[D(a)]„,™ ^ + (5„,™ (17) 

{1 71 Tfl 
' , is the Kronecker delta with ti, m G Z and the a„'s 

0, n ^ m 

arc the coefficients of the Fourier scries expansion of Q{u), 

oo 

Q{u)= ane^K (18) 



n— — OO 
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Thus, Eq. gives a stability criterion for the OHS mode by the condition 



1 - 2sin^ /C(K^)Vao det D(0) 



< 1, stable mode 
> 1, unstable mode 



(19) 



In this case, the Fourier coefBcients of Hill's matrix D(0) are given by the 
relations [Copson, 1935] 



ao 



2 -5 + 4^2 + 6 ^ ' '"="2 



TTJ 



drnr'^q" 



ICiK^)J \^2(7V + 1) 

• 2/ TTJ 



/C2(k2)(1 -qSn) \2{N+1) 



(20) 
(21) 



where q = e^'^^ , K. = JC{k^) and £ = £{k^) are the elliptic integrals of the first 
and second kind respectively and K.' = /C(k'^) with k'^ = 1 — k^. 

In the evaluation of D(0) in Eq. (|19|l we have used 121 terms in the Fourier 
series expansion of sn^(M, k^) (that is, 121 x 121 Hill's determinants of D(0)). 
Thus, we determine with accuracy 10"'^ the = k| values at which the Floquet 
exponent a in (|16|) becomes zero and the Zj{u) in H14|) become unbounded. We 
thus find that the first variation Zj{u) to become unbounded as increases is 
j = and the energy values at which this happens (see Eq. ((TT)) ') are 
listed in Table for j3 = 1.04. The Zj{u) variation with j = has Floquet 
exponent a equal to zero for every € [0, that is z n+i ju) corresponds to 
variations along the orbit. 



N 


5 


7 


9 


11 


13 


15 


Ec 


6.4932 


3.0087 


2.2078 


1.8596 


1.6669 


1.5452 



Table 1: The energies Ec, for f3 = 1.04, at which the z n-i (u) in H14|l becomes 
unbounded for some odd values of A^. 

Next we vary N and calculate the destabilization energy per particle for 
P = 1.04 at which the OHS nonlinear mode ^ becomes unstable. Plotting the 
results in Fig. ^ we see that decreases for large N with a simple power-law, 
as 1/iV. 



2.1.3 Analytical study of an OPM solution 

We now tu rn to the properties of another SPO of the FPU Hamiltonian Q, 
studied in [Budinskv fc Bountis. 1983t IPoggi fc Ruffo. 1997j and [Cafarella et 
al., 2003]. In particular, imposing the periodic boundary conditions 

XN+iit) ^ xiit), yt (22) 
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100 



500 



N 

Figure 1: The solid curve corresponds to the energy per particle , for (3 = 1.04, 
of the first destabilization of the OHS nonlinear mode © of the FPU system 
(0} obtained by the numerical evaluation of the Hill's determinant in H19|l . while 
the dashed line corresponds to the function oc Note that both axes are 
logarithmic. 



with N even, we analyze the stability properties of the out-of-phasc mode 
(0PM) defined by 

x,{t) = ~x,+i{t) = x{t), j = l,...,N. (23) 

In this case the equations of motion ||7J) reduce also to a single differential 
equation 

'i{t) = -4x{t) - 16/3x3 (i) (24) 

describing the anharmonic oscillations of all particles of the initial lattice. The 
solution of Eq. (|24|l can again be written as an elliptic cn-function 

x{t) ^Ccii{Xt,K^) (25) 

with 

The energy per particle of the nonlinear 0PM H23|l is given by 
in this case. 

We study the linear stability of the 0PM (|23|1 following a similar analysis 
to the one performed in the case of the OHS mode of Sec. 12.1.21 In this case, 
the corresponding variational equations have the form 

jjj = (1 + 12f3x^){y,-i - 2y, + v,+i), j^l,...,N (28) 



where yi = un+i- After the appropriate diagonaUzation, the above equations 
are transformed to a set of TV uncoupled independent Lame equations, which 
take the form 



+ {I + Ak^ -6K^sn'^{u,K^))sm^(^jzj{u) = 0, j = 1, . . . , TV (29) 



after changing to the new time variable u ~ At. Primes denotes again differen- 
tiation with respect to u. 

As in Sec. 12.1.21 Eq. H19|) gives a stability criterion for the nonlinear mode 
^ with 



ao 



-5 + An^ + 6^W[^ 



/C(k2) 



(30) 
(31) 



Proceeding in the same way as with the OHS mode, we find that the first vari- 
ation Zj{u) in Eq. (|29() that becomes unbounded (for = 1 and TV even) is ?' = 
^ — 1, in accordance with Budinskv fc Bountis. 1983| | , |Poggi fc Ruffo. 1997j , 
Cafarella et al. 2004| and that this occurs at the energy values E^, listed in 



Table 121 for (3=1. The Zj{u) variation with j = corresponds to variations 
along the orbit. 



TV 


4 


6 


8 


10 


12 


14 


Ec 


4.4953 


0.9069 


0.5314 


0.3843 


0.3051 


0.2532 



Table 2: The energies i?c, for /3 = 1, at which the zn__t^{u) in (|29|l becomes 
unbounded for the first time for some even values of TV. 

Taking now many values of N (even) and computing the energy per particle 
jf- for /3 = 1 at which the 0PM 1)23(1 first becomes unstable, we plot the results 
in Fig. 121 and find that it also decreases following a power-law of the form 
oc 1/TV2. 



2.2 The BEC model 

2.2.1 Analytical expressions for SPOs 

The Hamiltonian of the Bose-Einstein Condensate (BEC) model studied in this 
paper is given by 

N N N 

^ = 2 + '^^^ + I + ~ i Y.^P^P^+^ + m,+i) (32) 

i=i i=i i=i 



10 



100 



500 



N 

Figure 2: The solid curve corresponds to the energy per particle j^, for /3 = 1, 
of the first destabilization of the nonlinear 0PM of the FPU system (0} 
obtained by the numerical evaluation of the Hill's determinant, while the dashed 
line corresponds to the function oc Note that both axes are logarithmic. 



where 7 and e are real constants, which we take here to be 7 = e = 1. The 
Hamiltonian (|32|l possesses a second integral of the motion given by 

N 

P = J:(p" + i') (33) 

and therefore chaotic behavior can only occur for > 3. 

Imposing periodic boundary conditions to the BEG Hamiltonian (|32|l 

qN+i{t) = gi(i)and 

PN+i{t) = Pi{t), yt (34) 
we analyze the stability properties of the in-phase-mode (IPM) 
qj{t) EE q{t), 

Pjit) EE p{t) Vj = 1, . . . , TV, TV e N and ^ > 2 (35) 
and of the out-of-phasc mode (0PM) 
qj{t) = -qj+i{t) = q{t), 

Pj (t) = -Pj+iit) EE p{t) Vj = 1, . . . , iV, with N only even (36) 

and determine the energy range < i? < Ec{N) over which these two SPOs are 
linearly stable. 

In both cases, the corresponding equations of motion 

P3 = + J(p,2 + 92)q^.-l(g^._l+g^.+l)), Vj = l,...,7V (37) 
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give for the IPM solution 



k=V+l:{f + q^)p-^P, P=-(q+l:if +f)q-(q] (38) 



2 

and for the 0PM 

^ = p+^ip^ + f)p + ^p, -(^q+^{p^ + f)q + £qy (39) 

From Eq. (|33|l we note that the second integral becomes for both SPOs 

F^N{f+f) (40) 

yielding for the IPM solution 

^=(^-^+PjP^ P^-fl-+i^)^ (41) 



and for the 0PM 



2 / \ 2 



The above equations imply for both SPOs that their solutions are simple 
trigonometric functions 

q{t) = Ci cos{ujt) + C2 sin(wt) ^ (43) 
p{t) = — Ci sin(ct;t) + C2 cos{u}t) 

with w = 1 - e + ^ for the IPM and cj = 1 + e + 2£ for the 0PM with Ci and 
C2 real, arbitrary constants, where F = NA and A = Cf + C|. 
The energy per particle for these two orbits is then given by 

— = ^—^A + ^A^ (IPM) and — = + ^A^ (0PM). (44) 

N 2 8^' N 2 8^ ' ^ ^ 

Such SPOs have also been st udied in the case o f the integrable so-called 
dinier problem by other authors [Aubrv et al.. 199fij |. who were interested in 
comparing the classical with the quantum properties of the BEG Hamiltonian 
(El. 

2.2.2 Stability analysis of the SPOs 

Setting now 

qj ~ qj ^3 ' 

P] = Vj = 1,...,A^ (45) 
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in the equations of motion H37|l and keeping up to linear terms in xj and i/j we 
get the corresponding variational equations for both SPOs 



ij = (^-0 Vj-i + L+Vj - Vj+i + Kxj, 

Vj = " ((^l) ^^-1 + ^-^J - (l) ^^+1 + , = 1, . . . , TV (46) 

where xq = XN,yo = un and x^v+i — xi,yN+i ~ yi (periodic boundary condi- 
tions) and 



K 



L_ = 



1 



7(2Ccos(2wi) + Bsin(2wt)), 



1 



9, 



7(7j^ = 1 + A'j B7Cos(2a;t) + Cjsm{2Ljt) 



where B = C2 — and C = C1C2 are real constants. 

Unfortunately, it is not as easy to uncouple this linear system of differential 
equations (|46|l . as it was in the FPU case, in order to study analytically the 
linear stability of these two SPOs. We can, however, compute with arbitrarily 
accuracy for every given N the eigenvalues of the corresponding monodromy 
matrix of the IPM and 0PM solutions of the BEG Hamiltonian 

Thus, in the case of the 0PM (|36|l we computed for some even values of N 
the energy thresholds Ec at which this SPO becomes unstable (see Table 



N 


2 


4 


6 


8 


10 


12 


Ec 


5.0000 


4.5000 


3.1875 


2.4289 


1.9554 


1.6346 



Table 3: The energy Ec of the BEG Hamiltonian ^ at which the OPM ^ 
first becomes unstable for some even values of N. 

Plotting the results in Fig. 13 we observe again that Ec/N decreases with N 
following a power-law cx as in the case of the OPM of the FPU model. 

On the other hand, the IPM orbit (|35|l was found to remain stable for all the 
values of N and E we studied (up to = 54 and E w 10^). 

3 Destabilization of SPOs and Globally Chaotic 
Dynamics 

Let us now study the chaotic behavior in the neighborhood of our unstable 
SPOs, starting with the well-known method of the evaluation of the spectrum 
of Lyapunov Exponents (LEs) of a Hamiltonian dynamical system, Li, i = 
1, . . . , 2N where Li = imax > L2 > ■ ■ ■ > L2n- The LEs measure the rate 
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Figure 3: The solid curve corresponds to the energy per particle of the first 
destabilization of the 0PM of the BEG Hamiltonian (|22Jl obtained by the 
numerical evaluation of the eigenvalues of the monodromy matrix of Eq. (|46|) , 
while the dashed line corresponds to the function oc Note that both axes 
are logarithmic. 

of exponential divergence of initially nearby orbits in the phase space of the 
dynamical system as time approaches infinity. In Hamiltonian systems, the LEs 
come in pairs of opposite sign, so their sum vanishes, X]i=^i Li = and two 
of them are always equal to zero corresponding to deviations along the orbit 
under consideration. If at least one of them (the largest one) Li = Lmax > 0, 
the orbit is chaotic, i.e. almost all nearby orbits diverge exponentially in time, 
while if Lmax ~ the orbit is stable (linear divergence of initially nearby orbits) . 
Benettin et ai, [1980a, b] studied in detail the problem of the computation of 
all LEs and proposed an efficient algorithm for their numerical computation, 
which we use here. 

In particular, Li = Li{x{t)) for a given orbit x{t) is computed as the limit 
for t ^ oo of the quantities 



where Wi(0) and Wi{t), i = 1, . . . , 2N are deviation vectors from the given orbit 
x(t), at times t ~ and t > respectively. The time evolution of Wi is given by 
solving the so-called variational equations, i.e. the linearized equations about 
the orbit. Generally, for almost all choices of initial deviations Wi{0), the limit 
i oo of Eq. H49|l gives the same Li. 

In practice, of course, since the exponential growth of Wi{t) occurs for short 
time intervals, one stops the evolution of Wilt) after some time Ti, records the 
computed K?p_^ , orthogonormalizes the vectors Wi (t) and repeats the calculation 
for the next time interval T2, etc. obtaining finally Li as an average over many 



t \\m{0)r 

= lim Kj 



(48) 



(49) 
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Tj, j = 1,2, . . . ,n given by 



1 " 



(50) 



Next, we varied the values of the energy E keeping N fixed and studied 
the behavior of the Lyapunov exponents, using as initial conditions the OPMs 
(|23|) and of the FPU and BEC Hamiltonians respectively. First, we find 
that the values of the maximum Lyapunov exponent Li increase by two distinct 
power-law behaviors {Li cx , i? > 0) as is clearly seen in Fig.^lfor the 0PM 
of the FPU system. The result for the Li for the po wer-law behavior show n by 
solid line in Fig. lHis in agreement with the results in Rochester et oL. 1979j and 



Benettin. 1984j . where they obtain B = 0.5 and B = 2/3 for low dim ensional 



systems and differs slightly from the one obtained in Livi et al.. 1986j for the 
higher-dimensional case of = 80. We also find the same power-law behaviors 
with similar exponents B for the other positive Lyapunov exponents as well. 
For example, for the L2 we obtain L2 oc E^-^'^ and L2 oc E'^-^^ and L3 cx E^-^^ 
and La cx with the transition occurring at _B = 28.21 for all of them. 

FPU OPM N=10 



B,=0.63 
B =0.45 




50.0 62.5 75.0 



Figure 4: The two distinct power-law behaviors in the evolution of the maximum 
Lyapunov exponent Li as the energy grows for the OPM l)2;-{(l of the FPU 
Hamiltonian I^J for = 10. A similar picture is obtained for the L2 and L3 
also, with similar exponents and the transition occurring at the same energy 
value (see text). Note that both axes are logarithmic. 



Turning now to the full Lyapunov spectrum in Fig. [Slwe see that, for fixed 
iV, as the energy is increased (and more eigenvalues of the monodromy matrix 
exit the unit circle) the Lyapunov spectrum tends to fall on a smooth curve 
for the OPM orbits of FPU and BEC, see Fig. Ela), (b), as well as for the 
OHS mode with periodic boundary conditions (Fig. E^c)). Observe that in Fig. 
[SJ^c) we have plotted the Lyapunov spectrum of both the OPM (|23|l of the FPU 
Hamiltonian and of the OHS mode © for iV = 16 and periodic boundary 
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conditions at the energy E = 6.82 where both of them are destabihzed and 
their distance in phase space is such that they are far away from each other. We 
clearly see that the two Lyapunov spectra are almost identical suggesting that 
their chaotic regions arc somehow "connected" , as orbits starting initially in the 
vicinity of one of these SPOs visit often in the course of time the chaotic region 
of the other one. In Fig. [SJ^d) we have plotted the positive Lyapunov exponents 
spectra of three neighboring orbits of the OHS mode for N = 15 dof in three 
different energies and observe that the curves are qualitatively the same as in 
Fig. Ela) and (c). 



FPUOPMN=10 BECOPMN-10 




Figure 5: (a) The spectrum of the positive Lyapunov exponents for fixed = 10 
of the 0PM (123) of the FPU Hamiltonian Q as the energy grows, (b) Also for 
TV = 10, the 0PM ^ of the BEG Hamiltonian ^ yields a similar picture 
as the energy is increased, (c) The Lyapunov spectrum of the 0PM (|2;-{|) of 
the FPU Hamiltonian ijU for A'' = 16 and the OHS mode © of the same 
Hamiltonian and N, for periodic boundary conditions practically coincide at 
E = 6.82 where both of them are destabilized, (d) The Lyapunov spectrum 
of the FPU OHS mode © with fixed boundary conditions for iV = 15 as the 
energy grows presents as shape which is qualitatively similar to what was found 
for the SPOs of panel (c). 

Finally, in Fig. Elwc have plotted the eigenvalues of the monodromy matrix 
of several SPOs and have observed the following: For the 0PM of the FPU 
Hamiltonian the eigenvalues exit from —1 (as the orbit destabilizes via period- 
doubling) and continue to move away from the unit circle, as E increases further, 
see Fig. Eta). By contrast, the eigenvalues of the 0PM of the BEG Hamiltonian 
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exit from +1 by a symmetry breaking bifurcation and for very large E tend to 
return again to +1, see Fig. I3b). This does not represent, however, a return 
to globally regular motion around this SPO, as the Lyapunov exponents in its 
neighborhood remain far from zero. Finally, in Fig. El^c), we show an example 
of the fact that the eigenvalues of the IPM orbit of the BEC system, remain all 
on the unit circle, no matter how high the value of the energy is. Here, = 6, 
but a similar picture occurs for all the other values of N we have studied up to 
TV = 54 and E « 10^ 
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Figure 6: The eigenvalues , j = 1, . . . , 2iV (a) of the 0PM ^ of the FPU 
Hamiltonian ^ with N = 10 dof. (b) The eigenvalues of the 0PM ^ of the 
BEC Hamiltonian H32() with the same number of dof TV. (c) The eigenvalues of 
the 1PM (EH) of the BEC Hamiltonian with = 6 dof. 



4 Using SALI to Estimate the "Size" of Islands 
of Regular Motion 

In this section, we estimate the "size" of islands of regular moti on around sta- 
ble SPOs using the Smaller Alignment Index fSALI) metho d |Skokos. 20011 
ISkokos et ai. 200.Sat ISkokos et aL. 200.3bt ISkokos et al. 2004|. to distinguish 
between regular and chaotic orbits in the FPU and BEC Hamiltonians. The 
computation of the SALI has proved to be a very efficient method in revealing 
rapidly and with certainty the regular vs. chaotic nature of orbits, as it exhibits 
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a completely different behavior for the two cases: It fluctuates around non-zero 
values for regular orbits, while it converges exponentially to zero for chaotic 
orbits. The behavior of the SALI for regular motion was studied and explained 
in detail by |Skokos et ai. 2003bj , while a more analytical s tudy of the behavio r 
of the index in the case of chaotic motion can be found in ISkokos et ai. 2004j . 

As a first step, let us verify in the case of chaotic orbits of our A^-degree of 
freedom systems, the validity of SALI's dependence on the two largest Lyapunov 
expo nents Li = Lma.y an d L2 proposed and numerically checked for N = 2 and 
3, in [Skokos et aZ. . 20041 

SALI(t) (X e-(^i-^^)*. (51) 

This expression is very important as it implies that chaotic behavior can be 
decided by the exponential decay of this parameter, rather than the often ques- 
tionable convergence of Lyapunov exponents to a positive value. 

To check the validity of (|51|1 let us take as an example the OHS mode © of 
the FPU Hamiltonian Q using fixed boundary conditions, with A^ = 15 dof and 
f3 = 1.04, at the energy E = 21.6 and calculate the Lyapunov exponents, as well 
as the corresponding SALI evolution. Plotting SALI as a function of time t (in 
linear scale) together with its analytical formula (|51|1 in Fig. Cfa), we see indeed 
an excellent agreement. Increasing further the energy to the value E = 26.6878, 
it is in fact possible to verify expression H51|) . even in the case where the two 
largest Lyapunov exponents are nearly equal, as Fig. {7\h) evidently shows! 



FPU OHS N=15 



FPU OHS N=15 





Figure 7: (a) The time evolution of the SALI (solid curve) and of the Eq. H51|l 
(dashed line) at E = 21.6 of the OHS mode ® of the FPU Hamiltonian ^ 
with fixed boundary conditions, (b) Similar plot to panel (a) but for the larger 
energy E = 26.6878 for which the two largest Lyapunov exponents Li and L2 
are almost equal while all the other positive ones are very close to zero. In both 
panels the agreement between the data (solid curve) and the derived function 
of Eq. (|51|l (dashed line) is remarkably good. Note that the horizontal axes in 
both panels are linear. 



Exploiting now the different behavior of SALI for regular and chaotic orbits, 
we estimate approximately the "size" of regions of regular motion (or, "islands" 
of stability) in phase space, by computing SALI at points further and further 
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away from a stable periodic orbit checking whether the orbits are still on a 
torus (SALI> 10"^) or have entered a chaotic "sea" (SALR 10"®) up to the 
integration time t = 4000. The initial conditions are chosen perturbing all the 
positions of the stable SPO by the same quantity dq and all the canonically 
conjugate momenta by the same dp while keeping always constant the integral 
F, given by Eq. H33|) . in the case of the BEC Hamiltonian H32|) and the energy E 
in the case of the FPU Hamiltonian Q . In this way, we are able to estimate the 
approximate "magnitude" of the islands of stability for the 0PM of Hamiltonian 
(0} and for the IPM and 0PM of Hamiltonian H32|) varying the energy E and 
the number of dof N. 
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Figure 8: (a) "Size" of the islands of stability of the 0PM ^ of the BEC 
Hamiltonian (|32|l for = 4, 6 and 8 dof and SPOs constant energy E = 1 
before the first destabilization (see TableOI) . (b) "Size" of the islands of stability 
of the same Hamiltonian and SPO as in (a) for iV = 6 dof and three different 
energies of the SPO before the first destabilization (see Table IS)), (c) "Size" of 
the islands of stability of the same Hamiltonian as in (a) of the IPM (|35|1 for 
N = 6 dof and four different energies of the SPO. (d) "Size" of the islands of 
stability of the same Hamiltonian as in (a) of the IPM (ES) for f = ^. Here 
E corresponds to the energy of the IPM. 



In the case of the 0PM solutions of both Hamiltonians, as the number of dof 
N increases, for fixed energy E, the islands of stability eventually shrink to zero 
and the SPOs destabilize. For example, this is seen in Fig. |HIa) for the islands 
around the 0PM (|36(l of the BEC Hamiltonian (|32|l . A surprising behavior, 
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however, is observed for the same SPO, if we keep A'' fixed and increase the 
energy: Instead of diminishing, as expected from the FPU and other examples, 
the island of stability actually grows, as shown in Fig.|SJb), for the case of TV = 6 
dof. In fact, it remains of considerable size until the SPO is destabilized for the 
first time, through period-doubling bifurcation at E ^ 3.1875, whereupon the 
island ceases to exist! 

But what happens to the island of stability around the IPM solution of the 
BEG Hamiltonian H32|l . which does not become unstable for all values of N 
and E we studied? Does it shrink to zero at sufficiently large E or N7 From 
Fig-EJc) we see that for a fixed value of TV, the size of this island also increases 
as the energy increases. In fact, this SPO has large islands about it even if the 
energy is increased, keeping the ratio E/N constant (see Fig. ISJd)). This was 
actually found to be true for considerably larger E and N values than shown in 
this figure. 



Lyapunov Spectra and the Thermodynamic 
Limit 



Finally, choosing again as initial conditions the unstable OPMs of both Hamil- 
tonians, we determine some important statistical properties of the dynamics 
in the so-called thermodynamic limit of E and TV growing indefinitely, while 
keeping E/N constant. In particular, we compute the spectrum of the Lya- 
punov exponents of the FPU and BEG systems starting at the OPM solutions 
H23|l and H3t)|l for energies where these orbits are unstable. We thus find that 
the Lyapunov exponents are well approximated by smooth curves of the form 
Li « Lie""'/^, for both systems, with a w 2.76, a w 3.33 respectively and 
i = l,2,...,K{N) where K{N) « 3TV/4 (see Fig.©. 
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Figure 9: (a) Positive Lyapunov exponents spectrum of the OPM 1)23(1 of the 
FPU Hamiltonian I^J for fixed ^ = |- (b) Positive Lyapunov exponents spec- 
trum of the OPM (EH) of the BEG Hamiltonian ^ for fixed f = |. In both 
panels i runs from 1 to A^. 



Specifically, in the case of the OPM of the FPU Hamiltonian (g)) for 
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fixed energy density ^ = | we find 

L,(iV)«Li(7V)e-2-76^^ (52) 

wliile, in the case of the 0PM (|36|) of the BEG Hamiltonian (|32|l for fixed energy 
density ^ = f , a similar behavior is observed, 

L,(iV) «ii(7V)e-3-33^. (53) 

These exponential formulas were found to hold quite well, up to z = K{N) k, 
3-/V/4. For the remaining exponents, the spectrum is seen to obey different decay 
laws, which are not easy to determine. As this appears to be a subtle matter, 
however, we prefer to postpone it for a future publication. 

The functions H52() and (|53() . provide in fact, invariants of the dynamics, in 
the sense that, in the thermodynamic limit, we can use them to evaluate the 
average of the positive Lyapunov exponents (i.e. the Kolmogorov-Sinai entropy 
per particle) for each system and find that it is a constant characterized by the 
value of the maximum Lyapunov exponent Li and the exponent a appearing in 
them. 

In Fig. 1101 we compute t he well-known Kolmogorov-Sinai entropy husiN) 
Pesin. 197aTHilborn. 199i solid curves), which is defined as the sum of the 
A'^ — 1 positive Lyapunov exponents, 

hKs{N) = ^*(^)' ^'(^) > 0. (54) 

1=1 

In this way, we find, for both Hamiltonians, that h^siN) is an extensive thermo- 
dynamic quantity as it is clearly seen to grow linearly with TV (/iks(-/V) oc N), 
demonstrating that in their chaotic regions the FPU and BEG Hamiltonians 
behave as ergodic systems of statistical mechanics. 

Finally, using Eqs. (|52|l and (|53|) . as if they were valid for alH = 1, . . . , A — 1, 
we approximate the sum of the positive Lyapunov exponents Li, and calculate 
the hKs{N) entropy from Eq. H54() as 

hKsiN) (X Ln,,^ L,„ax « 0.095, (FPU 0PM) (55) 

— 1 + e N 

and 

/iKs(iV) « i,„ax L,„ax « 0.34, (BEG 0PM). (56) 

— 1 + e N 

In Fig. ^1 we have plotted Eqs. (|55|l and (|56|l with dashed curves (adjusting 
the proportionality constants appropriately) and obtain nearly straight lines 
with the same slope as the data computed by the numerical evaluation of the 
/iifs(7V), from Eq. ifK^ . 
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FPU 0PM E/N^3/4 BEC 0PM E/N=3/2 




Figure 10: (a) The hKs{N) entropy of the 0PM lO of the FPU Hamiltonian 
(0J for fixed ^ ~ j (sohd curve) and the approximated formula H55|l (dashed 
curve), (b) The hKsiN) entropy of the 0PM ^ of the BEC Hamiltonian 
for fixed ^ = | (solid curve) and the approximated formula (dashed 
curve) . 

6 Conclusions 

In this paper we have investigated the connection between local and global dy- 
namics of two N dof Hamiltonian systems, describing ID nonlinear lattices with 
different origins known as the FPU and BEC systems. We focused on solutions 
located in the neighborhood of simple periodic orbits and showed that as the 
energy increases beyond the destabilization threshold, all positive Lyapunov ex- 
ponents increase monotonically with two distinct power-law dependencies on 
the energy. We also computed the destabilization energy per particle of the 
OHS mode and of the two 0PM orbits and found that it decays with a simple 
power-law of the form E^/N oc iV~", a = 1 or 2. One notable exception is the 
IPM orbit of the BEC Hamiltonian, which is found to be stable for any energy 
and number of dof we considered! 

Furthermore, we found that as we increase the energy E of both Hamilto- 
nians for fixed N, the SPOs behave in very different ways: In the 0PM case of 
the FPU Hamiltonian the eigenvalues of the monodromy matrix always move 
away from —1 with increasing energy, while in the 0PM of the BEC Hamilto- 
nian these eigenvalues return to +1, for very high energies. Nevertheless, this 
behavior does not represent a return of the system to a globally regular motion 
around the SPO as one might have thought. It simply reflects a local prop- 
erty of the SPO, as the orbits in its neighborhood still have positive Lyapunov 
exponents which are far from zero. 

We have also been able to estimate the "size" of the islands of stability 
around the SPOs with the help of a recently introduced very efficient indicator 
called SALI and have seen them to shrink, as expected, when increasing N for 
fixed E in the case of the 0PM of the FPU Hamiltonian. Of course, when 
we continue increasing the energy, keeping the number of dof fixed, the OPMs 
destabilize and the islands of stability are destroyed. Unexpectedly, however, in 
the case of the 0PM of the BEC Hamiltonian these islands were seen to grow 
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in size if we keep N fixed and increase E up to the destabilization threshold, 
a result that clearly requires further investigation. The peculiarity of our BEC 
Hamiltonian is even more vividly manifested in the fact that the islands of 
stability about its IPM orbit never vanish, remaining actually of significant 
size, even for very high values of E and N (keeping E/N fixed). 

Starting always near unstable SPOs, we also calculated the Lyapunov spectra 
characterizing chaotic dynamics in the regions of the OHS and 0PM solutions 
of both Hamiltonians. Keeping N fixed, we found energy values where these 
spectra were practically the same near SPOs which are far apart in phase space, 
indicating that the chaotic regions around these SPOs are visited ergodically 
by each other's orbits. Finally, using an exponential law accurately describing 
these spectra, we were able to show, for both Hamiltonians, that the associated 
Kolmogorov-Sinai entropies per particle increase linearly with TV in the ther- 
modynamic limit of E —^ oo and N oo and fixed E/N and, therefore, behave 
as extensive quantities of statistical mechanics. 

Our results suggest, however, that, even in that limit, there may well exist 
Hamiltonian systems with significantly sized islands of stability around stable 
SPOs (like the IPM of the BEC Hamiltonian), which must be excluded from a 
rigorous statistical description. It is possible, of course, that these islands are 
too small in comparison with the extent of the chaotic domain on a constant 
energy surface and their "measure" may indeed go to zero in the thermodynamic 
limit. But as long as E and N are finite, they will still be there, precluding the 
global definition of probability densities, ensemble averages and the validity of 
the ergodic hypothesis over all phase space. Clearly, therefore, their study is 
of great interest and their properties worth pursuing in Hamiltonian systems of 
interest to physical applications. 
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8 Figure Captions 

1. The solid curve corresponds to the energy per particle for /3 = 1.04, of 
the first destabilization of the OHS nonlinear mode ^ of the FPU system 
Q obtained by the numerical evaluation of the Hill's determinant in (|19() , 
while the dashed line corresponds to the function cx Note that both 
axes are logarithmic. 

2. The solid curve corresponds to the energy per particle for /3 = 1, of 
the first destabilization of the nonlinear 0PM (|23|l of the FPU system Q 
obtained by the numerical evaluation of the Hill's determinant, while the 
dashed line corresponds to the function oc Note that both axes are 
logarithmic. 

3. The solid curve corresponds to the energy per particle jf- of the first 
destabihzation of the 0PM of the BEC Hamiltonian obtained 
by the numerical evaluation of the eigenvalues of the monodromy matrix 
of Eq. H46|) , while the dashed line corresponds to the function (x . Note 
that both axes are logarithmic. 
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4. The two distinct power-law behaviors in the evohition of the maximum 
Lyapunov exponent Li as the energy grows for the 0PM (|23() of the FPU 
Hamiltonian (0} for TV = 10. A similar picture is obtained for the L2 and 
L3 also, with similar exponents and the transition occurring at the same 
energy value (see text). Note that both axes are logarithmic. 

5. (a) The spectrum of the positive Lyapunov exponents for fixed iV = 10 of 
the 0PM of the FPU Hamiltonian igj) as the energy grows, (b) Also 
for N ^ 10, the 0PM ^ of the BEG Hamiltonian ^ yields a similar 
picture as the energy is increased, (c) The Lyapunov spectrum of the 
0PM (HSl) of the FPU Hamiltonian gj) for TV = 16 and the OHS mode 
© of the same Hamiltonian and TV, for periodic boundary conditions 
practically coincide at E = 6.82 where both of them are destabilized, (d) 
The Lyapunov spectrum of the FPU OHS mode ^ with fixed boundary 
conditions for A'^ = 15 as the energy grows presents as shape which is 
qualitatively similar to what was found for the SPOs of panel (c). 

6. The eigenvalues A^, j = 1, . . . , 2A (a) of the 0PM ^ of the FPU Hamil- 
tonian (gj) with A^ = 10 dof. (b) The eigenvalues of the 0PM ^ of the 
BEG Hamiltonian (|32|l with the same number of dof A^. (c) The eigenval- 
ues of the IPM (123 of the BEG Hamihonian ^ with A^ = 6 dof. 

7. (a) The time evolution of the SALI (solid curve) and of the Eq. H51|l 
(dashed hne) at E = 21.6 of the OHS mode ® of the FPU Hamiltonian 
g} with fixed boundary conditions, (b) Similar plot to panel (a) but 
for the larger energy E — 26.6878 for which the two largest Lyapunov 
exponents Li and L2 arc almost equal while all the other positive ones are 
very close to zero. In both panels the agreement between the data (solid 
curve) and the derived function of Eq. ^31^ (dashed line) is remarkably 
good. Note that the horizontal axes in both panels are linear. 

8. (a) "Size" of the islands of stability of the 0PM ^ of the BEG Hamilto- 
nian 132|l for A = 4, 6 and 8 dof and SPOs constant energy E = 1 before 
the first destabilization (see Table ^ . (b) "Size" of the islands of stabil- 
ity of the same Hamiltonian and SPO as in (a) for A' = 6 dof and three 
different energies of the SPO before the first destabilization (see Table 13). 
(c) "Size" of the islands of stability of the same Hamiltonian as in (a) of 
the IPM ^ ioT N = 6 dof and four different energies of the SPO. (d) 
"Size" of the islands of stability of the same Hamiltonian as in (a) of the 
IPM ^ ior § = ^. Here E corresponds to the energy of the IPM. 

9. (a) Positive Lyapunov exponents spectrum of the 0PM (|23(l of the FPU 
Hamiltonian gj for fixed ^ = |. (b) Positive Lyapunov exponents spec- 
trum of the 0PM lEH) of the BEG Hamiltonian ^ for fixed f = |. In 
both panels i runs from 1 to A^. 

10. (a) The hKs{N) entropy of the 0PM ^ of the FPU Hamiltonian Q 
for fixed ^ = f (solid curve) and the approximated formula H55|l (dashed 
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curve), (b) The hxsiN) entropy of the 0PM ^ of the BEG Hamiltonian 
(|32(l for fixed ^ = f (sohd curve) and the approximated formula H56|l 
(dashed curve). 
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